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ABSTRACT 

This study follows the development and initial employment of piezo- 
resistive accelerometers as ship motion sensors. The spectral results 
are discussed as they relate to concurrent measurements of wind fluctua- 
tion and temperature variance data. Both analog and digital Fourier 
analysis were performed. Results include co-spectra, quad-spectra, and 
phase and coherence plots. Biasing of mean parameters by pitching effects 
are also discussed. 
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I. INTRODUCTION 



The purpose of this study was to develop measurement procedures for 
shipboard motion and to analyze resulting data for information about 
platform motion effects on turbulence measurements. The shipboard motions 
of interest included the large scale motion associated with the waves, 
intermediate scales associated with heavy equipment vibration, and high 
frequency motions associated with vibrations of the instrument mast. 

The results of this study are necessary in evaluations of different 
approaches for making turbulence measurements from a ship. In the past, 
turbulence measurements have been primarily made from land stations or 
stable overwater platforms. Open-ocean turbulence data are scarce be- 
cause of the difficulties in describing turbulence properties from ship- 
board measurements. 

This study is a part of a larger investigation being conducted by 
an interdisciplinary group at the Naval Postgraduate School on optical 
wave (laser) propagation in the marine boundary layer. The primary 
objective of the group is to relate optical propagation to the measure- 
able properties of atmospheric turbulence. The area of concentration is 
to define the properties and effects unique to the marine environment. 
Since these effects depend upon turbulence properties related to the 
wind, temperature, humidity and wave conditions, a subgroup was formed 
to conduct micrometeorological observational experiments. 

The purpose of the boundary layer investigation is to gather the 
data necessary to define the state of atmospheric turbulence during the 
propagation experiments. For this it is necessary to make both tur- 
bulence and profile measurements from a ship. During these experiments 

the ship is located along the propagation path. 
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The scarcity of actual data that can be used to verify existing 
boundary layer theory requires further empirical examination. While 
present boundary layer theory is primarily based on overland results, 
there is reason to believe that the turbulent regime of the marine 
boundary layer may differ significantly from that found in land-based 
studies, due to water wave effects. 

This study is directly concerned with measurement problems encountered 
from unstable platforms, specifically ships. The phases in the study 
were the development of the measurement system, the measurement of motion 
during a turbulence experiment, and analysis of the results. 
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II. TURBULENCE MEASUREMENT CONSIDERATIONS 



Parameters used in boundary layer descriptions include the wind 
shear stress parameter (U*) and the heat flux parameter (T*) . The first 
parameter is a direct derivative of the mean flux of horizontal momentum 
in the vertical direction (-u'w 1 ) which is. related to the vertical wind 
shear (du/dz). The second is also related to the vertical temperature 
gradient (de/dz). The emphasis in this study is primarily on the wind 
shear stress parameter (U^) . 

The wind shear stress parameter appears in the following expressions 
for the wind profile and wind gradient: 



U = |* [Ln(z/z o )J 

and 

dU _ U* 
dz kz 



( 1 ) 

( 2 ) 



where U is the mean horizontal wind at height (z) over the surface of 

characteristic roughness length (z ). The constant of proportionality 

o 

(k) is the Von Karman constant. 

Generally these expressions and other similar expressions have 
been formulated or evaluated with actual measured fluxes. This has been 
done almost exclusively from land experiments or from stable marine plat- 
forms. Such comparisons cannot be made from shipboard systems due to 
the unacceptable error induced by tilt angles as small as one degree 
to the vertical axis when measuring vertical fluxes. 

However, on the basis of the Kolomogorov hypothesis for the spectral 
distribution of velocity variance and on the basis of the turbulent 
kinetic energy balance expression, one can obtain estimates of the 



7 



dissipation of turbulent kinetic energy (e) and the dissipation of 
temperature variance ( e g) which can then be related to the fluxes. 

Dissipation estimates are part of the Kolomogorov 1 s hypothesis, i.e. 
there exists, in wave number space, a separate region between the fre- 
quencies which generate and dissipate turbulent kinetic energy. In 
this region, referred to as the inertial subrange, neither generation 
nor dissipation occurs. Turbulent kinetic ‘ energy is transferred to 
higher frequency disturbances until it is dissipated by viscous forces. 

From this hypothesis, a -5/3 power law for the inertial subrange of 
wind fluctuation spectra was suggested (Figure 1). The equation has the 
form 



+ (k) - C 



2/3 k -5/3 



(3) 



where ^(k) is the power spectrum of wind fluctuations; e is dissipation 
of turbulent kinetic energy; k is the wave number of the disturbance 
within the inertial subrange, and C is a constant. 

There is another relation for estimating . It is based on the 
wave number (k^) which corresponds to the upper frequency limit of the 
inertial subrange as illustrated in Figure 1. The relation between and 
the wave number corresponding to the upper limit of inertial subrange is: 
k U 3 . 

~~ = C*-) 174 (4) 

2 7T e 



Both relations should be fully exploited for best results, and both 
require accurate description of spectral properties. Further discussion 
is available in Johnston (1974). 

The dissipation rate ( e) can also be related to stability parameters 
and the stress parameter (U^) by the equation: 



e 




(-) - -] 
m '‘I/ L J 



(5) 
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Figure 1. Turbulence Regime in Wavenumber Space. 
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where ^^ (— ) is an empirically derived stability function dependent on 
the air-sea temperature difference. L is the MONIN-OBUKOV stability 
length. 

Similar expressions exist which relate the dissipation of temperature 

variance ( e ) to larger scale thermal properties of the boundary layer, 

0 

but are not discussed here. 

The key in an approach in which flux estimates are obtained from 
dissipation estimates is the ability to obtain accurate descriptions of 
the spectra associated with wind fluctuations (u f ) , and temperature 
fluctuations (t T ). With good spectral data reliable estimates of e and 
can be made. Initial experiments aboard the R/V Acania yielded veloc- 
ity spectra which appeared contaminated due to ship-related motion and 
thus demonstrated the need for this study. Figure 2 illustrates the 
results of unpublished studies by Kasales (1973) aboard R/V Acania. The 
recurrent spikes and improper slope of produced spectra required the 
development of a sensor system to check the origin of such spectral 
deviation. 
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Figure 2. Kasales Wind Fluctuation Spectra. 
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III. DEVELOPMENT OF MEASUREMENT SYSTEM 



A. SUMMARY OF DEVELOPMENT 

The experiment, which commenced on 13 September 1973 and lasted ten 
days, was the final phase of a continuous program of observational re- 
search and sensor development which began early in 1973. The plan to 
examine ship motion was initiated in July 1973. Sensor availability 
with only two months of lead time proved to be a limiting factor in 
the experiment. ‘Other meteorological sensors such^ as quartz thermo- 
meters were in procurement processes for six months or more. Some 
had been used in previous experiments and on the earlier cruise to San 
Nicholas Island. Other sensors and systems were being used on ship- 
board experiments for the first time. 

A field test of a piezoresis tive accelerometer was performed in 
early August after preliminary sensor evaluation. The reason for the 
test was to check the sensor 1 s very low frequency response and also to 
estimate its signal to noise level. For the test, one piezoresis tive 
accelerometer was taken aboard the R/V Acania for an afternoon cruise in 
Monterey Bay. The test was hindered somewhat due to spray associated 
with pitching in large swell conditions. However, sufficient data were 
recorded to verify a satisfactory low frequency response to the wave 
field and to verify that there was sufficient signal level. The particu- 
lar accelerometer used in this test had less than half the rated signal 
level of the sensor ultimately used, and also had a higher resonant 
frequency. 
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The selection of the desired sensor, and the construction of the 
conditioning unit were begun after completion of the shipboard feasi- 
bility trial. Procurement was expedited yet all development tasks were 
completed only slightly more than a week before the experiment. 

Calibration was limited due to time available and the necessity for 
syttem check-out. Calibration information from tests by the manufacturer 
on each of the six instruments was available and were depended on heavily. 
In particular, response changes with frequency had to be accepted. How- 
ever, the manufacturer’s signal level calibrations could not be used 
due to a change in the excitation voltage supplied to the accelerometer 
by a modified conditioning unit. Calibration of signal levels were 
completed prior to installation aboard the R/V Acania. 

B. SELECTION CRITERIA 

Several criteria were defined as being important in the final selec- 
tion of the instrument to be used as the motion sensor. Factors which 
later proved to be most important were signal output level, cost, and 
availability as well as range of frequency response. Both gyroscopes 
and accelerometers were included in the initial consideration. Gyro- 
scopic systems would have been most valuable in interpreting low fre- 
quency motion. They also would have provided inherent stable orientations 
of the axes. However, they were either too large or too expensive for 
use in this experiment. A typical cost would have been in excess of 
thirty thousand dollars for a triaxial prepackaged system, and a lease 
would have been about five thousand dollars for the necessary time frame. 

After tests on the sensitivity of the instrument on board the R/V 
Acania, (described previously), the piezoresis tive accelerometer was 
chosen. Its advantages were its very small size, weight , and linear 
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response over a broad frequency range. The piezoresistive transducer 
consists of a seismic mass and four semiconductor crystals which support 
the mass, two in tension and two in compression, Figure 3. The resist- 
ance of the crystalline material changes with a change in stress on the 
crystal. In a bridge circuit, the resistance change produces a signal 
voltage. Any acceleration of the mass along its axis of support by the 
crystals will cause a change in crystal stress. The lack of dependence 
on crystal support except along a single axis provides directional sen- 
sitivity with only a small transverse sensitivity. The resulting signal 
voltage in a bridge circuit is directly proportional to the acceleration, 
and is sensitive even to very low frequency or steady state accelerations 
such as gravity. 

The thermal effects on the sensitivity of the accelerometer were not 
as critical as first envisioned. Thermal effects must be considered 
whenever a resistance type sensor is being used. Tests indicated a de- 
viation of less than three percent from the measurement of the manufac- 
turer’s rated (75°F) temperature over a range from 0 to 160°F. During 
the experiment the maximum temperature variation was nine Fahrenheit 
degrees which occurred over several data periods. A thermal zero shift 
due to temperature changes was noticeable in analog traces extending over 
periods of hours. This was associated with diurnal temperature variations. 

Gravity contamination proved to be the most serious problem, as was 
expected. Depending on the axis of orientation the change in the com- 
ponent if gravitational acceleration due to the axial tilting would ex- 
ceed the value of the translational acceleration. In particular, a beam 
sea produces roll of the R/V Acania often in excess of ten degrees, hence 
the component of the gravitational acceleration can exceed the actual 
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Mass 



Figure 3. Model Diagram for a Piezoresis tive Accelerometer. 

t 
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lateral acceleration. However this effect is expected to occur only 
within a narrow low frequency spectral band. Because of the gravity 
effect, pitch is generally a more reliable axis for low frequency ac- 
celerometer data. This is fortunate since the ship was headed into 
the wind and sea during measurements in order to obtain better sampling 
of the wind field. Thus, the motion due to pitch was greater than that 
due to roll. Gravity contamination limited quantitative analysis of wave 
and ship interaction frequencies for this study. This could probably be 
overcone by different mounting arrangement, such as mentioned by Miyake 
( 1970 ). 

Integration of the accelerometer data was not a serious problem. 
Analog integration was performed and the results were excellent for low 
frequency study. Such integration is a form of low pass filtering. This 
would seem to nullify the advantage of the broader frequency spectrum, 
however, this is not the case. Assuming a relatively small number of 
sinusoidal Fourier components, which when integrated keep the same fre- 
quency but shift in phase, it does not matter if the accelerations are 
integrated or not for spectral analysis studies. 

The piezoresis tive accelerometers were selected in part because of 
their availability. Fortunately those available from the supplier 
covered the expected frequency range and had a high signal to accelera- 
tion ratio. The problem bandwidth versus signal level was minimized. 

In this regard, the minimum bandwidth necessary should be such that the 
signal level will sufficiently exceed noise. 

It is characteristic of the device that the resonant frequency (F^) 
determines the bandwidth. Two easy estimates of acceptable bandwidth 
are the four percent error point of one-fifth of resonant frequency, 
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and the half power (-3 db) point of one-half resonant frequency. Graph- 
ical plots of frequency response normalized by resonant frequency (rela- 
tive frequency response) can be found in Appendix A. 

Another important consideration was the susceptibility of the sensors 
to shock damage due to dropping or metal to metal contact. In fact, one 
of the six sensors used in this study failed in spite of special shock 
blocking included by the manufacturer, and careful handling by the ex- 
perimenter. 

Without regard to the long-term zero shift previously mentioned, the 
reproducibility of measurements was better than one percent in calibration 
runs and approached three digit accuracy for four sensors out of six. 

This calibration accuracy exceeded that of all instruments used in the 
experiment with the exception of the quartz thermometers. Post experi- 
mental calibration was not successful due to recording and conditioning 
gain problems during the experiment. 

C. CONDITIONING OF SIGNALS 

The accelerometers required additional components before voltage sig- 
nals could be recorded. Two primary considerations were that the signal 
had to be amplified from a few millivolts up to the order of one volt, 
and that the Wheatstone bridge circuit had to be completed. Although a 
system was available from the manufacturer, construction of a bridge 
system for six accelerometers was accomplished at the Naval Postgraduate 
School. 

The Wheatstone bridge in this system was located in two sepatate 
parts and joined by the cables. The side of the bridge corresponding 
to the accelerometer changes resistance (in the manner previously de- 
scribed) in response to the excitation voltage allowing a varying current 
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to flow through one-half of the bridge circuit. The other half of the 
bridge circuit consistes of a matched pair of 1.7 K ohm resistors. 

Additional bridge components consisted of a 10 K ohm current limiting 
resistor, a balance potentiometer for nulling the steady state current 
and the excitation power supply. The power supply vas solid state con- 
trolled and had a direct current output of five volts. 

The manufacturer^ rated excitation voltage for the accelerometers 
was twice that supplied by the conditioning unit which was used. This 
resulted in a signal level degradation by a factor of nearly three from 
that specified for the sensor. Instead of signal levels of over 30 
millivolts per g (gravity force) unamplified, voltages were of the order 
of 10 to 12 millivolts per g. Bridge signals were amplified by a solid 
state millivolt amplifier with resistors matched for a hundredfold gain, 
but limited by a variable resistor used as an output attenuator. The 
signal voltage was that which existed between this resistor and ground. 

It was amplified by gains of from 30 to 60 before recording. 

A diagram of the accelerometer-cab le- conditioning unit is shown as 
Figure 4. Unfortunately, the conditioning unit was inadvertently de- 
calibrated during a tape recorder change. This also was complicated 
by recorder calibration problems. 

D. MOUNTING CONSIDERATIONS 

The accelerometer triads were mounted during the experiment at the 
base of the mast 2.54 m above the waterline and at the highest level of 
meteorological instruments approximately 14.42 m, Figure 5. The lower 
mounting was envisioned primarily as a pitch and roll mounting with 
the hope that it would also measure any specific hull vibrations trans- 
ferred to the mast through its mountings. The upper mounting was designed 
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Figure 4. Circuit Diagram for the Measurement System. 
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Figure 5. Vane Assembly for Upper Level, Mast Mounting at Lower Level. 
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to describe as exactly as possible the vibrations to which the wind 
fluctuations probe was subjected, with the secondary objective of showing 
by comparison with the lower sensor, any mast structure deviations from 
solid body motion. Few deviations due to mast motion were expected for 
low frequencies because of guy wires attached to the mast at two levels. 
The lattice structure of the mast was expected to contribute some unique 
high frequency vibrations not found in the spectra of the lower sensor. 

Special considerations were made to insure proper mounting of the 
accelerometers. To give the greatest reproduction of motion the mounting 
system had to be as rigid as possible. In this regard the upper triad 
of sensors were mounted in an aluminum casing and directly connected to 
the wind fluctuation probe assembly by tapping screws into the base block 
of the probe assembly. This was done to maximize the correlation between 
the vibrations sensed by the wind fluctuation probe and those sensed by 
the accelerometer group. Aluminum was chosen as the material for the 
container since both the size of the casing and its weight could affect 
the wind resistance and balance of the vane assembly which aligned the 
probe with the mean wind. To minimize these effects a second casing of 
the same size with added internal weight was added to the opposite side 
of the vane probe assembly. The vane with both casings attached was at 
the highest level on the mast. Accelerations at that level were compared 
with the wind fluctuation spectra at that level as well as with acceler- 
ations measured at deck level. Figure 6. 

The lower sensor package was constructed of dressed one- inch lumber 
and was eight inches square and nine inches deep. It was waterproofed 
with a jam lid, silicone rubber sealant and waterproof canvas tape. 
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MOUNTING 
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Mounting was on angle iron about four inches off the deck surface to 
prevent pooling of water, and additionally rigidized with nylon strapping 
tape. The lower triad was aligned with the ship’s centerline. 

Internally, both mountings were similar with emphasis on rigidity 
and thermal insulation. The purpose of the latter was to reduce the 
effect temperature changes associated with solar radiation, relative 
wind shifts which would change the ventilation of the casing surface, 
and possibly turbulent temperature fluctuations. One ‘procedure used 
to reduce temperature changes was to paint both casings white. Addition- 
ally each set of three sensors was mounted on an acrylic block with one 
sensor on each of three mutually perpendicular faces. 

Acrylic plastic is solid and is therefore good for vibration trans- 
mission. It is also a poor heat conductor which provided extra thermal 
protection over metallic blocks. The surface of the acrylic can be 
buffed very smooth and suffers less deformation damage than wood surfaces. 
The buffing provided a good contact and minimal warping of the accelero- 
meter face. Unfortunately the acrylic proved to be too brittle to 
allow tapping holes for mounting screws so the accelerometers were mounted 
on double sided tape as a compromise rather than permanently affix them 
with epoxy glue. This measure proved satisfactory because of the small 
sensor weight along with the additional support from styrofoam pellets 
which filled the cavities providing the final thermal insulation. The 
acrylic block was then glued to the outer casing, and the casings were 
sealed. 

Cabling was attached permanently with waterproofed connections to 
each accelerometer. The cable was a t riaxial-grounded shield type. 

Cables were 200 feet in length. The cables from upper and lower sensor 
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packages were strung from the foredeck of the R/V Acania to the dry 
laboratory where they were connected to the conditioning unit. Cables 
were color coded with respect to sensors and conditioning amplifiers 
during calibration and were not changed during the experiment. 

Initial shipboard checkout during the very calm conditions of early 
morning indicated that one sensor had failed during installation. There- 
fore the deck mounted triad of sensors was reoriented to exclude the 
fore and aft accelerations in favor of pitch and roll accelerations. 
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IV. THE EXPERIMENT 



A. SUMMARY OF THE INSTRUMENT SUITE 

The observational experiment was conducted aboard the R/V Acania, 
Figure (7), operated by the Department of Oceanography. The ship is over 
48 meters in length, has 6.5 meter beam and has a draft of 2.7 meters. 

The complete instrument suite consisted of groups of meteorological 
instruments mounted at several levels on a 50-foot mast on the foredeck 
of the ship and a short mast at the bow. Additionally, a sea-surface 
temperature probe was trailed from a boom amidships. 

The optical portion of the experiment included a gyro-stabilized 
helium neon laser mounted in the wet laboratory. The laser receiving 
equipment was located on the northwest tip of San Nicolas Island. 

The mean parameter components of the meteorological instrument suite 
included cup anemometers for mean wind measurements, hygrosensors for 
humidity measurements, and quartz thermometers for mean temperature read- 
ings. Additional information on these sensors has been provided by 
Cavanaugh (1974). 

Fluctuations in temperature and velocity were measured with hot-wire 
type sensors used with different overheat characteristics. Further dis- 
cussion of these is provided by Johnston (1974). Photographs of instru- 
ment arrays appear in Appendix B. 

B. RECORDING OF DATA 

Conditioned signals were recorded on a Sangamo FM 14-channel tape 
recorder. Difficulties with each of three tape recorders on board the 
R/V Acania during the experiment reduced the amount of usable data to 
small portions of each data tape and, additionally, limited the longest 
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Figure 7 . 



R/V Acania with meteorological mast 
(Stem mast not shown) 
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data period available to 20 minutes. However, sufficient periods of data 
were available from a variety of recording periods with different mean 
conditions to meet most experimental objectives. 

A difficulty encountered after the data were gathered was that the 
recorder used could not reproduce the data on all channels at all avail- 
able speeds since playback amplifiers did not have the speed circuitry 
required. In addition an intermittent capstan drive oscillation caused 
loss of large sections of data due to large spurious noise signals which 
blanketed the recorded data. Finally, it was discovered that tape 
recorder channels had different gains associated with them, and had not 
been marked for proper replacement so that exact calibration of signal was 
impossible . 

These conditioning and recording problems greatly reduced the quanti- 
tative value of the remaining data. A primary problem during the experi- 
ment proved to be equipment malfunction, which may occur during shipboard 
use of sophisticated electronic equipment. 

C. THE OBSERVATIONAL PHASE 

The transit from Monterey Bay to San Nicolas Island was the beginning 
of the data collection segment of the experiment. During the transit the 
instrument signals were monitored on an oscilloscope. Winds from the 
stern blew exhaust gases over the sensor arrays during this phase, and 
precluded measurement of the turbulent regime. This effect and super- 
structure wind interference greatly limit the allowable directions of 
relative wind for representative turbulence measurements while underway. 

A strip chart record was made of the output of the five functioning 
accelerometers during the transit, prior to the arrival at San Nicolas 
Island . 
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The ship anchored on the night of 19 September 1973, in 10 fathoms of 
water about 1600 yards from the shoreline of San Nicolas Island. The 
anchorage was on the northwest side of the island and allowed unlimited 
fetch to the winds which were from the northwest, Figure (8). The kelp 
moderated the smaller wind waves, and water depth probably increased the 
swell height. 

On the morning of 20 September 1973, data collection from all systems 
was begun. Unfortunately, the tape recorder used for the accelerometer 
data began to act erratically and attempts to locate the problem were 
unsuccessful beyond isolation of the problem to a caixd in the power sup- 
ply. The problem was later determined to be a cracked component. 

Fortunately channels were available on the recorder used to record 
wind and temperature fluctuation data and they were \used. Coincident 
recording enabled examination time correlation and phase information 
between these two sets of data. If possible, this should be done on 
future studies. 

D. SUMMARY OF MEASUREMENTS 

Acceleration measurements were made with five sensors (one failed 
very early in the experiment) mounted in groups of three and two at the 
two levels. All three upper level sensors and the lower roll and pitch 
sensors lasted through the five days of constant use and rough sea con- 
ditions. Due to large direct current offset voltages with respect to 
signal output much data was lost due to the fact that the tape recorder 
was preset with insufficient voltage limits. Fortunately, the usable 
data were those which were recorded simultaneously with velocity and 
temperature data (when available). 

The time frames and parameters available as well as descriptions of 
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Figure 8. North End, San Nicolas Island 
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mean conditions (from Cavanaugh, 1974) appear in Appendix C. Only a small 
portion of the available data was analyzed due to long analysis procedures 
necessary to proceed from the original analog signals to the final digi- 
tal computer output. 

E. ANALYSIS PROCEDURES 

Analysis procedures consisted of a post-cruise strip playback of all 
analog tapes for preliminary analyses of the data, digitalization . of 
selected time periods, digital processing of the data from seven- to nine- 
track tape, and Fourier analysis to yield spectra. 

In the first step, analog signals were played back at high speeds and 
the output were displayed on an eight-channel strip— chart recorder run at 
a very slow speed. This yielded an overview of the data available for 
analysis. One segment of such a recording is included as Figure (9) 
where data include wind fluctuation (u ! ), temperature fluctuation (t ! ), 
and accelerometer data (A). Coincident mean data were recorded by hand 
at 10-minute intervals. This provided 10-minute averages of winds, with 
median values for the other variables. 

Time segments with usable data were restricted due to the recording 
problems previously mentioned. The length of continuous data periods was 
limited to about 20 minutes and precluded sequential studies of overlap- 
ping time frames. 

Data were reproduced with the same tape recorder used to record the 
data, but only after some repairs had been effected. Prior to conversion 
a gain of 50 was applied to the input signals, and one channel of the 
accelerometer data was integrated. 

The integration was done in order to relate low frequency ship related 
motion in a relative velocity to fluctuation measured by the velocity 
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Figure 9. Fluctuation variable signals. From left, 

3(u f ), 2(t ! ), 2(A), and integrated accelerometer. 
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fluctuation sensors . 



During the digital conversion, velocity signals were low-pass filtered 
at 400 Hertz in order to remove spurious spikes. Direct current offsets 
were also nulled. It is thought that an integrated signal for one time 
frame (run 3) was influenced by an improper null adjustment. 

One of the problems associated with multiple processing of data is the 
addition of noise. High frequency noise is added due to non-linear elec- 
trical responses. Power supply noises and their harmonies are also cumu- 
lative. Spurious spikes can result from connectors and these multiply as 
the number of equipment interfaces increase. All such noise signals con- 
taminate the final spectra. A direct input to an analog fast Fourier 
transform device (spectrum analyzer) would lessen many of these problems. 

Analog to digital conversion was performed on the Naval Postgraduate 
School SDS 9300 analog/digital hybrid computer system. These procedures 
have been detailed by Jones (1971) and Johnston (1974) and yield a seven- 
track digital tape. This must be converted to a nine- track tape for use 
in the IBM 360 system digital computer. Several minutes of computer 
time are required per conversion. Data from this experiment produced 13 
digital tapes. The seven- to nine-track tape conversion is described in 
the W. R. Church Computer Center, Technical Note No. 0211-08. 

The nine-track digital tape was processed with a version of the 
University of British Columbia Fourier Transform program which has been 
adapted for use at the Naval Postgraduate School. It is used to compute 
complex Fourier coefficients which are used to describe spectra evergy 
distribution. Final output requires another program which takes the 
tabular coefficients as input and converts them to power and cross- 
spectra which are then output in final form on a Cal-Comp plotter. 
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Johnson (1974) has described the procedures used on this experiment to 
obtain the final spectra. 
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V. RESULTS 



A primary outcome of this experiment was that the developmental 
system functioned properly throughout the experiment, giving data which 
could be compared to the fluctuation parameters. The results from the 
observational experiment include strip chart recordings, analog variance 
spectra, and digital variance spectra, cospectra, quad-spectra as well 
as phase and coherence spectra. These various results were examined 
for correlations between the accelerometer data and the wind fluctuation 
data. 

A. STRIP CHART RESULTS 

Strip chart results were useful for the interpretations. This fact 
is illustrated by the strip chart sections in Figure 9 which show accel- 
erometer outputs, at low chart speed, in comparison with the wind fluct- 
uations (u T ) and temperature fluctuations (t ! ) signals. Figure 10 
illustrates the relation between signals fron the five accelerometers. 

An examination of the same signal reproduced at higher chart speeds 
shows that high-frequency components were different between sensors, 
Figures 11 and 12. 

Other strip chart recordings suggest significant correlations between 
the single low frequency wind fluctuation signal available and the inte- 
grated vertical accelerometer signal. This relationship is most evident 
for larger motion oscillations, and less significant for small motions. 
Figure 13 illustrates this phenomenon for a section of the strip chart. 
The correlation is further evident between the vertical accelerometer 
data and wind fluctuation data both of which had been high-pass filtered 
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Figure 11. Accelerometer signal (inertial frequencies) lOmm/sec. 
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Figure 13. Wind fluctuation data compared with integrated upper vertical acceleration data. 



at ten Hz. Although the signal was high-pass filtered there is still 
evidence of some energy due to low frequency ship's motion. This is only 
evident for very large signal levels. Figure 14 shows this low frequency 
coupling. 

The sections of strip-chart recording which were included in this 
section were those which revealed the low frequency motion effects on 
the wind fluctuation sensor. There were only two strip chart records 
with the proper gain and speed relationships to make the evidence visible. 
Unfortunately, low frequency correlation analysis of spectra could not 
be accomplished due to a combination of the digitization rate and block 
size limitations in the computer code. Spectral results are described 
in the following section. 

B. ANALOG SPECTRAL RESULTS 

Analog spectra were obtained from an x-y plotter connected to a 
SPECTRAL DYNAMICS MODEL SD-301C Real Time Analyzer and SD-309 Ensemble 
Averager, available during the cruise. The x-y plotter was used to make 
spectral plots i n nearly real time. Many such plots made on 20 Sep- 

tember 1973. Unfortunately, for the available spectra, the system was 
not properly calibrated so that scales are not consistent. However, 
decade values could be marked so frequency can be estimated. Variance 
spectral plots of the velocity (u 1 ) obtained by this method are presented 
in Figures 15, 16, 17, and 18. 

Accelerometer spectra from the upper triad appear in Figure 19. It 
is interesting to note the coincidence of spectral spikes among the three 
sensors. Especially significant is the fact that the predominant spike 
is many decibels larger than the secondary spikes. The frequency associ- 
ated with the predominant spike is 54 Hertz. A similar plot of spectra 
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Figure 14. Filtered wind data compared with upper vertical acceleration data. 




Figure 15. Analog spectrum of u' , level (1) . 
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Figure 17. Analog Spectrum of u' , level (3). 
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Figure 18. Analog spectrum of u' , level (4). 
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Figure 19. Analog Spectra of Upper Triad Accelerometer Signals. 



for the lower two sensors does not have the striking resemblance of the 
upper triad plot. However, dissimilar spikes and differences in shape 
are evident in Figure 20. A tabular listing of significant spikes ac- 
cording to frequency is given in Table I below. 

TABLE I 

FREQUENCIES OF 500 HZ SPECTRUM 

264.5 120.0 30.0 

211.5 105.0 24.5 

173.5 54.0 16.5 

158.5 39.5 11.5 

Spectral differences are evident in combined spectral plots of the 
accelerometers and the u ? sensors. The predominance of spikes found 
in the accelerometer spectrum is barely discemable in the wind variance 
spectrum. The predicted decreasing spectrum of the inertial subrange 
with its associated -5/3 slope is clearly seen in Figure 21. 

C. DIGITAL SPECTRAL RESULTS 

Digital spectral results differed slightly in appearance from ana- 
log spectral results due to the different bandwidths in the two compu- 
tations. Digital spectra were obtained for three data periods in 
order to make visual comparisons between spikes in the wind variance 
spectra and accelerometer spectra. The accelerometer spectra have the 
same form for all these runs; the upper triad sensors with the analog 
integration of the vertical accelerometer as the fourth spectrum. Com- 
bined spectra are illustrated in Figure 22. 
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Figure 20. Analog Spectra of Pitch and Roll Accelerometer Signals. 
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Figure 21. Analog Spectral Shapes (not to same scale) . 



LOG , 0 SPECTRUM 
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Figure 22. Combined Digital Accelerometer Spectra. 
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In comparison with the analog spectra, there is less representation 
of detail in the digital spectra. Although the detail could be obtained 
by digital means, it would be at a significant cost in computer time 
and core storage. It may be further contraindicated since larger band- 
widths would tend to smooth any sharp vibration induced spike. Digital 
spectra for the data runs (2), (3), and (8) appear as Figures 23, 24, 
and 25, respectively. 

The latter run, Figures 25a - d, which contained the only unfiltered 
u T data, was selected for most detailed analysis. The objective was 
to see if functional relationships could be established. For this 
examination, cospectra, quad-spectra, phase, and coherence spectra were 
computed. These are interpreted below for phase relationships between 
the three accelerometer signals as well as the wind and temperature 
variance data. 

D. ADDITIONAL DIGITAL SPECTRA 

Three wind variance spectra and three accelerometer apectra appear 
in Figures 26 to 31. From the co-spectra, quad-spectra, phase and co- 
herence plots, conclusions can be drawn about the relation of the tur- 
bulent wind variance and the accelerometer data. It is noted that the 
integrated acceleration shows strong relationship with the wind fluctua- 
tion on basis of the co-spectra, quad-spectra, and coherence spectra for 
the lowest frequency band. Such similarity is notably absent in the 
region from 30 to 50 Hz in all four plots, concurring with the idea 
that dissipation rate (e) can be estimated at these wave numbers without 
much ship motion effect. Figures 32 to 35 illustrate this effect. 

Additionally the phase shows a remarkably constant phase, near zero, 
for frequencies above 60 Hz, while the coherence values rise almost 
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Figure 24. 
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Figure 25c. 
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Figure 25d. 
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Figure 27. 
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Figure 28. 
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Figure 29. 
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Figure 31. 
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Figure 32. 
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Figure 34. 



65 




U 3 LF.UEL 1 

INTEGRATED UPPER LF'.UEL UERT RCCEL 



Figure 35 



linearly with frequency beyond 60 Hz. This is interpreted as either 
noise intrusion into both spectra, or the effect of the wind energy 
spectrum decreasing (-5/3 law) until the motion effect gradually becomes 
more evident. 

Consideration of these same plots for the upper triad of accelero- 
meters firther supports a conclusion of isotropy at high frequencies, 
Figures 36 to 43. The similarity of the co-spectra and quad-spectra 
shapes as well as the phase and coherence plots for the lateral and 
vertical accelerometers show that a single instrument could probably 
suffice for all but the most exact description of the motion. The fore- 
aft motion and its relationship to vertical motion is more complex. The 
co-spectra plot shows a strong similarity with the previous co- 
spectral plot, but the quadrature spectrum show two distinct minima near 
80 and 200 Hertz which also show on the coherence plot. A possible 
explanation for this is vibrations due to rotating machinery aligned 
fore-aft. These vibrations would be sensed by both vertical and lateral 
accelerometers but not the one aligned with the shaft. 
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Figure 38. 
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Figure 40. 
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VI, CONCLUSIONS 



A, ABOUT THE MOTION MEASUREMENT SYSTEM 

The accelerometer system developed in this study was found to be 
sufficiently adapted for shipboard application. 

B. ABOUT TURBULENCE MEASUREMENTS FROM SHIPS 

The results of the experiment indicate that small scale turbulence 
measurements can be made from a ship. This is a qualified result due 
to the lack of data during low wind speeds. During the experiment mean 
winds ranged from about 10 to 20 knots with the latter being most fre- 
quent. These velocities are sufficient to give well developed turbulence 
and a fairly broad inertial subrange. It is apparent from Figures 45 
and 46 showing turbulent wind fluctuation spectra and accelerometer 
spectra that the inertial subrange is generally unaffected by motion 
between 20 and 50 Hertz. Thus the viscous dissipation parameter (c) 
can be estimated from this spectral region. Part of this insensitivity 
to motion may have been due to the differing nature of the wind fluctua- 
tion spectral form and the spectral form of the accelerometers. 

At frequencies above and below the inertial subrange the motion of 
the platform affected the turbulence measurements. Low frequency com- 
ponents in the strip chart segments indicate that there was certainly 
some fictitious wind induced as the sensor followed the ship f s pitching 
motion through the wind profile. 

The results by Kasales (1973) show a consistent energy peak in the 
vicinity of two Hertz. This was one of the aspects to be explored. It 
was at one time thought that this was due to a vane assembly oscillation. 
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However, this was not the case since two analog spectra from this study 
show similar energy peaks in the region of five Hertz. A significant 
fact was that the accelerometers had a minimum of energy in the same 
frequency band. Figure 44 is one of the spectra. This experimenter can 
only conclude that is a real property of atmospheric turbulence, but not 
a consistent one. 

The high frequency end of the spectrum was dominated by a single 
spike near 54 Hz, which is thought to have been the frequency of the 
shipboard power supply. This same contamination appears in u T spectra 
Figures 45 and 46. This noise, however, is not very significant since it 
occurs at a distinct frequency. 

Other high frequency noise affects the dissipation end of the inertial 
subrange and ultimately limits a vital cross check on the value of the 
stress parameter (U^) , as discussed in the background section. It is 
difficult to select the proper wave number (k Q ) for viscous drop-off 
with such noise in the spectra. Shock mounting of the sensors would 
decrease the high frequency contamination. 

Spectral analysis from both analog or digital imeans revealed that 
though the amplitude varies greatly for different members of the upper 
triad the frequencies remain the same. There were temporal variations 
in frequency of the spikes, but spatial consistency among the triad 
appeared to hold in general. With respect to spectral consistency between 
sensors accelerometers at each level appear to be sufficient to define 
the high frequency deviations at the dissipation end of the inertial 
subrange. 
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Figure 44. 
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Figure 45. 
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Figure 46. 



C. CONSIDERATIONS ON MEAN MEASUREMENTS FROM SHIPS 



The mean sensors also showed an effect due to the ship’s motion. The 
pitching of the vessel is the logarithmic wind profile and could produce 
the observed deviations. The greatest deviations occur at lower heights 
and would result in an apparent lowering of the sensor. Figure 47. The 
profile results described by Cavanaugh (1974) showed consistently high 
values of the drag coefficient, C^. It is possible that an empirical 
correction could be determined for this based on the measured shear and 
the magnitude of the vertical motions which cause this effect. 

In shipboard studies of these parameters Denman (1973) , did not ex- 
amine this factor. It is felt that the increase in ship motion during 
increasing winds may well cause over-estimates of shear, resulting in 
larger drag coefficients. 

D. SUGGESTIONS FOR FUTURE RESEARCH 

Specific change recommendations are to increase excitation voltage in 
the bridge circuit top to the rated level and to use single vertically 
stabilized accelerometers at each sensor level. 

A bottom mounted pressure sensor for wave field data would be use- 
ful in describing ship response to the wave field, and also the wave 
field effects on turbulence. Digitization should be performed to ex- 
amine low frequency information. 

An analog spectrum analyzer should be acquired and used to obtain im- 
mediate spectra. This would allow in situ evaluation of the data with 
more data runs. 



81 




Figure 47. Sensor height variation produces biased mean 
wind and wind shear values. 
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APPENDIX A 



SENSOR CHARACTERISTICS AND CALIBRATION RESULTS 



TYPE: Piezoresis tive accelerometer, 

MANUFACTURER: Endevco 

CHARACTERISTIC DATA: 

Excitation 
Range 

Shock limits 
Frequency response 
Weight 

Transverse sensitivity 
Warmup time 
Resistance (each arm) 

Ground 
Signal level 

Linearity and hysteresis 
Humidity protection 
FACTORY CALIBRATION: 21 August 1973 



Serial Color Signal Maximum 

Transverse 

Sensitivity 

BC 25 Red 34.4 mv/g .3% 

BC 50 White 33.5 1.9 

BC 65 Blue 36.8 3.0 

BC 28 Cromate 37.5 2.4 

BC 18 Grey 33.0 1.6 

BC 43 Gold 31.4 2.5 



Model 2265-20' half bridge type. 

10 volts DC 

±20 g 

+ 80 g, 2000 g pulsed maximum 
+ 5%, steady state to 200 Hz. 

6 grams 
Less than 5% 

1 minute 

Nominal 1800 ohnis 

Inner shield connected to cable shield 
Nominal 30 millivolts /g 
+ 2% over range 

Epoxy and silicone rubber sealed 



Zero 

Measured 

Output 


Resonant 

Frequency 


+13.5 mv 


1.11 KHz 


+14.3 


1.13 


+ 8.8 


1.12 


-19.7 


1.08 


+19.0 


1.08 


+ 9.5 


1.13 
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FREQUENCY RESPONSE: 





50 Hz 


100 Hz 


150 Hz 


200 Hz 


Red 


0 


0 


+1% 


+2.5% 


White 


0 


0 


+1 


+2 


Blue 


0 


0 


+1 


+2.5 


Chromate 


0 


0 


+1 


+2.5 


Grey 


0 


0 


+1 


+3 


Gold 


-.5 


0 


+1 


+2 



PRESAIL CALIBRATION (with amplifier): 

The presail calibration was limited to a signal output level measure- 
ment. The accelerometers were connected to 200 feet of cable and both 
were color-coded to a specific amplifier. The bridge was then balanced 
to +0.0001 volt on a recently calibrated Simpson digital voltmeter. 
Sensors were attached to a vane and positioned to measure the earth T s 
gravitational field. Measurements were also taken at ninety degree in- 
crements. These measurements were all completed at least twice. 









TABLE A-l 






AMP 


COLOR 






Volts 








LEVEL (H/L) 
AXIS 


"lg 


AV 


0g 


AV 


+lg 


6 


Red 

H.L. 

(Lateral) 


+0.027 

+0.028 


0.347 


-0.320 

-0.319 


0.378 


-0.698 

-0.697 


5 


White 


+0.205 




-0.188 




-0.605 




L.L. 


+0.208 


0.395 


-0.187 


0.416 


-0.603 




(Pitch) 


+0.211 




-0 . 186 




-0.602 


4 


Blue 

H.L. 

(Vertical) 


+0.063 

+0.063 


0.329 


-0.266 

-0.267 


0.370 


-0.636 

-0.635 
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TABLE A-l (con) 



AMP 


COLOR 




volts 








LEVEL (H/L) 
AXIS 


-lg 


AV 


Og 


AV 


+lg 


3 


Gray 

L.L. 

(failed) 


-0.030 

-0.030 


0.344 


-0.374 

-0.373 


0.394 


-0.769 

-0.768 


2 


Chromate 

H.L. 

(fore-aft) 


+0.475 

+0.478 


0.629 


-0.153 

-0.153 


0.646 


-0.799 

-0.798 


1 


Gold 

L.L. 

(roll) 


-0.196 

-0.192 


0.391 


+0.200 

+0.195 


0.387 


+0.582 

+0.586 



INSTALLATION CHECK-OUT: 

On completion of installation sensors were checked for operation 
during a calm night with R/V Acania moored pierside prior to getting 
underway. 

SHIP'S CLINOMETER READING 



List: 0.5 degree starboard 






Roll: 1 degree peak 


TABLE A- 2 




Sensor 


Mean Voltage 


Deviation 


Red 


0.243 


0.004 


White 


0.371 


0.001 


Blue 


0.240 


0.002 


Gray 


2.485 


Steady 


Chromate 


0.195 


0.002 


Gold 


0.258 


0.003 
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Figure A-l. Schematic of a piezoresistive accelerometer (from manufacturer). 




Figure A-2 . Phase shift as a function of damping (from manufacturer). 
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RELATIVE RESPONSE , Response - dB 




.05 .1 .2 .5 1.0 2.0 

Normalized Frequency 




5 

Frequency 

Figure A-3. Frequency response curves (from manufacturer). 
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figure A-5. 



Upper mounting after disassembly. 



Arrows indicate sensors. 
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Figure A-6 . Typical sensor. 



91 




92 






APPENDIX B. 



PHOTOGRAPHIC RECORD OF THE EXPERIMENT 




Figure B-l. Dry laboratory array. From left, spectrum analyzer, tape 
recorder, fluctuation conditioning panels, oscilloscope monitor. 
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Figure B-2 . Dry laboratory second view. 
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Figure B-3. Dry laboratory array. From left, accelerometer conditioning 

unit, filters, tape recorder. 
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Figure B-4. Dry laboratory array. Items included; 
monitor, humidity array monitor, wind speed counter 



temperature 
indicators . 



array 
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Figure B-5. LT W.E. Johnston checks probe assembly on stem mast. 
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Figure B-6. LT M.P. Cavanaugh climbs the main instrument mast to 

Connect cables. 
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Figure B-7. Highest level of mast, note accelerometer and blank casings 

on probe assembly. 
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Figure B-8. Lower mounting. 
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Figure B-9. 



Sea surface temperature probe. 
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Figure B-10. Close-up of wind speed counter indicator. 
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APPENDIX C 



SUMMARY OF MEASUREMENTS AND ANALYSIS 



DIGITIZED DATA PERIODS: 



TIME 




*FLUCTUATI0N VARIABLES 


**MEAN VARIABLES 


DIGITAL TAPES 


2000 Sept. 


20 


u' (4) , A, INTA 


u(3),t(4),T o ,q(3) 


2 


2323 Sept. 


20 


u' (4), A, INTA 


u(3),t(4),T o ,q(3) 


2 


0637 Sept. 


21 


u'(4), A, INTA 


u(3),t(4),T o ,q(3) 


2 



* A = Accelerometer data 

INTA = Integrated upper vertical accelerometer 

u 1 (3) = Wind fluctuation at three levels 



** u(2) = mean (10 minute) wind at two levels 
t(4) = mean temperature at four levels 



T q = sea surface temperature 
q(2) = humidity at two levels 



SPECTRAL ANALYSIS: 



TIME 




VARIABLES 


2000 


Sept. 20 


u'(4), A(3), 
INTA 


2323 Sept. 20 


u'(4), A(3), 
INTA 


0637 


Sept. 21 


u'(4), A(3), 
INTA 



REMARKS 

Run 2. Spectra of inertial subrange 
and accelerometers (Stable). 

Run 8. Spectra of inertial subrange 
and accelerometers by both analog and 
digital means. u| was the only very 
low frequency u 1 signal. Selected 
for co-spectra, quad-spectra, phase, 
and coherence studies. (Neutral) 

Run 3. Spectra of Inertial subrange 
and accelerometers. Integrated ac- 
celerometers suffered poor null 
adjust (Unstable). 
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